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ABSTRACT 

An  approximate  Boltzrnann  equation,  known  as  the  single 
relaxation  model  [k] ,   is  studied  here.  This  equation  is 
linearized  and  the  fundamental  solution  is  considered.  Fol- 
lowing H.  Grad,  [6],  the  solution,  asymptotic  in  small 
values  of  the  ratio  of  mean-free-path  to  distance  from  the 
origin,  is  sought.   It  can  be  shown  that  the  fundamental 
solution  itself  gives  the  asymptotic  description  of  the  flow 
field  past  an  object.  This  solution  gives  the  asymptotic 
description  when  the  distance  from  the  origin  is  much  greater 
than  either  the  mean-free-path  or  the  body  size.  This  is  true 
independently  of  the  Knudsen  ntmiber. 

It  is  shown  that  within  the  linearized  framework  used,  the 
Navler-Stokes  equations  give  a  consistent  asymptotic  description 
of  the  correct  kinetic  theory  flow  field.   On  the  other  hand, 
the  inviscid  Euler  equations  fail  in  the  neighborhood  of  the 
Mach  lines.   The  asymptotic  fundamental  solution  found  gives 
the  fine  structure  of  the  wake  uniformly  in  the  Mach  number. 
In  addition,  the  solution  gives  a  detailed  description  of  the 
neighborhood  of  the  Mach  line,  uniformly  in  Mach  number,  as  the 
latter  varies  from  supersonic  to  high  subsonic  values.   The 
uniform  solution  is  given  in  closed  form  in  terms  of  an 
Airy  function. 
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In  the  limit  of  zero  rr.ean-free-path,  the  fundamental 
solution  cf  the  invlscid  equations  is  recovered.   Con- 
vergence to  the  inviscld  solutions  is  not  uniform. 
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ON  THE  KINETIC  THEORY  OF  STEADY  GAS  FLOWS 

1.   Introduction. 

In  recent  years  there  has  been  an  Increasing  interest 
in  the  problems  of  kinetic  theory.  Much  of  this  renewed 
interest  can  be  traced  to  the  work  of  Wang  Chang  and 
Uhlenbeck  [7],  [14]^  and  later  Mott-Smith  [20].   Using 
approximate  Boltzmann  equations  and  exploiting  the  tools  of 
linearized  theory j,  they  were  able  to  deal  with  some  special 
but  fundamental  problems;,  such  as  sound  propagation,  Couette 
flow,  and  heat  flow„   More  recently  attention  has  been  paid 
to  kinetic  models.  By  a  kinetic  model  we  mean  an  equation 
which  is  intended  to  mimic  the  Boltzmann  equation  in  as  many 
ways  as  possible ^  and  yet  remain  analytically  simple.  A  model 
which  seems  to  have  first  appeared  in  the  literature  in  a  work 
of  Welander  [1]  and  simultaneously  in  a  paper  by  Bhatnager, 
Krook  and  Gross  [2],  has  proven  successful  in  a  n\jmber  of 
problems,  see  [2],  [4],  [5],  The  modeling  lies  in  the  replace- 
ment of  the  Boltzmann  collision  operator  by  a  term  proportional 
to  the  deviation  of  the  distribution  from  the  locally 
Maxwelllan  distribution  function 

The  constant  of  proportionality,  v,  is  the  collision  frequency 1 
the  remainder  of  the  notation  is  that  commonly  used  and  is 
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further  explained  in  Section  3.   This  equation  exhibits  the 
more  important  features  of  the  full  Boltzmann  equation. 
Briefly,  moments  of  the  collision  term  vanish  when  weighted 
by  the  collisional  Invariants,  so  that  the  conservation 
equations  follow,  furthermore  the  H  theorem  is  valid ^  for 
details  see  [8].   It  is  also  clear  from  the  model  that  a 
relaxation  theory  (the  space  homogeneous  case)  predicts  that 
all  moments  of  the  distribution  fiinctlon^  other  than  density, 
velocity,  temperature,  decay  with  the  same  period.  For  this 
reason,  the  equation  is  referred  to  as  the  single  relaxation 
model  [4].   Quite  recently  H.  Grad  [6]  considered  the  funda- 
mental solution  of  the  single  relaxation  model c   He  developed 
a  technique  for  gaining  an  asymptotic  form  of  the  fiindamental 
solution.  This  is  fully  discussed  in  Section  3- 

Before  going  on  further,  we  digress  a  moment  to  discuss 
what  is  meant  by  the  linearized  Boltzmann  equation.   As  is 
shown  In  [4],  our  equation  can  be  derived  from  the  Boltzmann 
equation  in*  which  squares  of  the  perturbed  distribution 
function  away  from  the  local  Maxwellian  are  neglected.  This 
equation  is  often  referred  to  as  the  "linearized"  Boltzmann 
equation.   Yet,  it  is  far  from  being  linearized,  since  moments 
of  the  distribution  function  occur  in  the  local  Maxwellian 
(a  quadratic  non-linearity  is  replaced  by  an  exponential  non- 
linearity).   We  shall  refer  to  this  linearization  as  the 
weak  linearization,  and  reserve  the  appellation  of  strong 
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linearization  for  a  kinetic  equation  which  is  linearized  away 
from  an  absolute  Maxwellian.  Aside  from  the  reference  cited 
above,  another  recipe  for  obtaining  the  single  relaxation 
model  (from  the  exact  non-linear  Boltzmann  equation  in  this 
case)  can  be  found  in  [3].  No  consideration  of  these  deriva- 
tions will  bt  found  in  this  workj  rather  if  one  wishes,  the 
merits  of  the  model  may  be  evaluated  in  light  of  the 
subsequent  theory  it  predicts. 

In  Section  J>,   a  procedure  which  resembles  the  Chapman- 
Enskog  method  [16],  is  applied  to  the  steady  strongly  linearized 
form  of  our  equation  (also  see  the  discussion  in  the  next 
section).  The  hierarchy  of  linearized  hydrodynamical  equations 
is  found  to  be:  the  linearized  Euler  equations  (Sec.  5)> 
the  Oseen  equations  (see  Section  5)*  a  system  formally 
equivalent  to  the  linearized  Biornett  equations,  and  higher 
moment  systems  (see  Section  6).   As  is  indicated  in  Section  5, 
the  computed  viscosity  coefficient  is  "correct"  ,  however,  the 
heat  conductivity  is  off  by  a  factor  of  J>/2.      It  may  seem 
rather  remarkable  that  the  Navier-Stokes  system  should  result 
after  so  much  has  been  done  to  destroy  the  character  of  the 
collision  operator.  The  implication  is  that  the  Naview-Stokes 


The  viscosity  (or  the  conductivity,  but  not  both)  can  be 
made  anything  one  wishes  by  the  appropriate  choice  of  v,  the 
"collision  frequency".   A  consistent  approach,  as  is  given 
in  [4]  indicates  the  choice  of  v.  This  leads  to  a  viscosity 
coefficient,  which  Chapman- Cowling  [l6],  call  the  first 
approximation. 
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system  is  relatively insensitive  to  the  collision  process, 
and  depends  on  only  those  gross  properties  which  have  been 
Included  in  the  single  relaxation  model. 

In  the  above  introduction,  we  have  anticipated  some  of 
the  results  which  will  be  found  later.   In  the  section  which 
follows,  in  addition  to  interpretation,  we  give  a  verbal 
pictiore  of  the  results  of  this  report.   A  detailed  study  of 
the  stjructure  of  certain  simple  flows  is  given  as  well  as 
application  to  more  complicated  flows.   No  attempt  at  proving 
statements  is  made;  however,  reference  to  later  discussion  is 
given  parenthetically. 


2„   Survey ♦ 

As  mentioned  in  the  preceeding  section,  a  kinetic  model 
will  be  used  in  place  of  the  full  Boltzmann  equation.   However, 
the  preliminary  remarks  made  in  this  section  apply  also  to  the 
full  Boltzmann  equation.   For  this  reason  we  start  the  dis- 
cussion in  generality  and  only  later  particularize. 

In  order  to  sharpen  the  succeeding  discussion,  let  us 
focus  attention  on  a  particular  problem.  We  consider  the  flow 
past  a  particle  source  located  at  the  origin.   The  source  can 
be  thought  of  as  a  point  in  space,  which  spews  out  and  absorbs 
particles.  To  further  fix  the  picture,  we  take  the  particle 
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source  such  that  it  conserves  the  number  of  particles,  i.e., 
it  is  not  a  mass  source.   Further,  we  endow  it  with  a  lift, 
drag,  and  heat  generation.   Higher  moments  such  as  stress, 
heat  flow,  etc.,  could  also  be  associated  with  the  particle 
source,  but  these  will  not  enter  the  discussion  which 
follows.   If  the  equation  governing  the  flow  is  taken  to  be 
linear,  the  proper  superposition  of  these  sources  can  be 
found  (in  principle)  to  give  the  flow  past  some  arbitrarily 
shaped  body.   For  our  strongly  linearized  equation,  this  is 
the  problem  of  finding  the  fundamental  solution  (see  Section 
3).  For  the  present,  it  is  not  necessary  to  be  this  restrictive. 

The  Boltzmann  equation  for  a  simple  gas  contains  only  one 
pertinent  length  scale,  the  mean  free  path,  L.  We  gain 
tractabillty  if  instead  of  seeking  a  general  solution,  we  look 
for  a  solution  which  is  asymptotic  in  either  small  l/R  or  small 
r/l.   R  is  the  distance  from  the  origin.   A  solution  correspond- 
ing to  the  latter  limit  describes  the  flow  field  within  a 
mean  free  path  of  the  origin,  where  we  have  essentially  "free 
molecular"  flow.   An  application  of  this  solution  in  the 
linear  case  would  be  the  superposition  of  sources  along  a  wall 
to  give  a  description  of  the  slip  flow  region  adjacent  to  the 
wall.   Incidentally,  in  this  regime  the  collision  term  of  the 
full  Boltzmann  equation,  which  contains  the  non-linearity,  is 
a  higher  order  term,  so  that  the  superposition  of  sources 
would  presumably  still  give  an  accurate  picture. 

-  10  - 


It  Is  the  former  limits  that  of  small  L/R,  with  which 
we  will  be  concerned  in  the  main  body  of  this  work.   A  solution 
corresponding  to  this  limit  depicts  the  flow  field  many  mean 
free  paths  from  the  origin.   Since  the  gas  many  mean  free  paths 
from  boundaries  has  had  the  opportunity  to  relax,  the  distri- 
bution function  will  be  nearly  Maxwellian.  We  then  woiild 
expect  the  deviation  from  the  Maxwellian  to  be  given  in  terms 
of  the  density,  temperature,  velocity,  and  their  derivatives. 
That  the  description  will  be  hydrodynamical  (i.e.,  in  terms 
of  density,  temperature,  and  velocity)  may  be  seen  in  another 
way.   With  proper  normalization  of  the  steady  Boltzmann 
equation,  the  Chapman-Enskog  expansion  is  an  expansion  in  the 
ratio  L/R,  supposedly  valid  for  small  L/R.  This  method,  as 
is  well  known,  leads  to  a  hierarchy  of  fluid  dynamical 
equations.   Although  the  language  used  is  fluid  mechanical, 
it  is  Important  to  note  that  the  above  remarks  refer  to  the 
entire  Knudsen  niwiber  range.  When  sufficiently  far  away  from 
a  body,  a  hydrodynamical  description  may  be  given,  independent 
of  the  rarity  of  the  gas. 

At  this  point  we  leave  the  general  case  to  place  restric- 
tions on  the  governing  equation  and  the  cases  considered.  We 
choose  the  steady  form  of  the  strongly  linearized  single 
relaxation  model  (see  equation  (3»10)).  The  constant  flow  at 
negative  infinity  is  governed  by  the  absolute  Maxwellian  about 
which  linearization  takes  place.  Linearization  places  a 
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rather  strong  restriction  on  a  solution.   We  know  from 
experience  that  linearization  Is  incorrect  within  shock 
waves  and  the  Prandtl  boundary  layer.  To  avoid  repetitious 
reference  to  the  linearization  requirement,  we  adopt  the 
point  of  view  that  we  are  dealing  with  a  gas  which  behaves 
in  accordance  with  the  linearized  equation.  Then  at  the 
end  of  the  discussion  we  can  make  the  a 11 -encompassing 
declaration  that  the  solution  is  not  valid  in  regions  of  the 
flow  where  quadratic  terms  are  not  negligible  with  respect 
to  linear  terms. 

Under  the  restrictions  of  the  preceeding  paragraph, 
the  problem  set  forth  is  equivalent  to  finding  the  fundamental 
solution  (see  Section  J>) .     H.  Grad  showed  that  to  the  first 
approximation  (in  L/r),  this  reduces  to  finding  the  fundamen- 
tal solution  of  the  linearized  Euler  equations.   In  anticipa- 
tion of  what  follows,  it  will  serve  our  purposes  if  we  pause 
to  consider  what  is  meant  by  first  approximation,  secor.d 
approximation,  etc.  The  first  approximation  comes  merely  from 
the  grouping  of  terms  of  0(L/r)  and  ignoring  the  higher  orders 
(see  Section  5);  the  second  approximation  consists  of  grouping 
terms  of  0(L/r)  and  0(L  /r  )  and  ignoring  the  rest  (see 
Section  5),  and  so  on.  This  procedure  results  in  a  hierarchy 
of  equations  in  density,  temperature  and  velocity.  As 
mentioned  in  the  introduction,  this  results  in  the  linearized 
Euler  equations,  the  Oseen  equations,  and  then  higher  moment 
systems. 
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Although  grouping  terms  in  this  manner  seems  computation- 
ally attractive,  there  is  no  compelling  reason  to  do  so. 
Rather  than  a  rule,  it  should  act  as  a  guide.   We  now 
elucidate  this  comment.   Solving  for  the  fundamental  solutions 
of  the  linearized  Euler  equations,  we  get  the  first  approxima- 
tion to  the  fundamental  solution  of  the  strongly  linearized 
single  relaxation  model.  This  leads,  in  the  case  of 
supersonic  flow,  to  singularities  along  the  Mach  cone,  see 
[6]  (also  see  Sec.  k).     Further,  H.  Grad  shows  in  reference 
[6],  that  there  is  another  singularity  which  is  along  the  flow 
direction  behind  the  source.  When  properly  interpreted,  this 
can  be  identified  with  the  wake  in  the  inviscld  limit  (Sec.  4). 
The  latter  singularity  exists  through  subsonic  velocities  (Sec. 
4).   In  the  neighborhood  of  these  singularities  higher  order 
terms  in  (l/r)  are  of  the  same  effective  magnitude  as  those 
terms  in  the  first  approximation  (Sec.  6).   For  this  reason 
the  grouping  according  to  0(L/R)  is  not  valid  throughout  the 
flow  field.   One  finds  that  the  terms  which  must  be  carried 
to  correctly  describe  the  "singular  regions"  (they  are  no 
longer  singular  with  the  additional  terms)  belong  to  the 
second  approximation,  i.e.,  they  come  from  the  grouping  which 
yields  the  Oseen  equations  (Sec.  5).   Furthermore,  the  terms 
added  by  higher  approximations  are  indeed  of  higher  order, 
and  in  particular  this  is  true  for  the  terms  which  come  from 
the  "Burnett"  equations  (Sec.  6).  This,  within  the  linearized 
kinetic  theory  framework  laid  down  here,  substantiates  the  use 

-  15  - 


of  the  Navler-Stokes  equations  in  fluid  mechanics.  Actually 
the  association  with  hydrodynamlcal  equations  need  never 
have  been  made.  In  the  sense  that  It  has  no  bearing  on  the 
analysis.   One  can  formally  carry  out  the  analysis  even  if 
in  complete  Ignorance  of  fluid  mechanics  (at  least  in  principle) 
It  should  again  be  noted,  that  although  we  have  been  using 
the  language  of  fluid  mechanics  (implying  small  Knudsen 
number),  the  results  and  analysis  apply  to  the  entire  spectrum 
of  Knudsen  number.  This  is  so  for  the  simple  reason  that  the 
fundamental  solution  is  applicable  to  bodies  of  any  size  and 
hence  independent  of  the  usually  defined  Knudsen  niomber.  There 
is  however,  if  we  choose  to  define  it,  a  Knudsen  nymber  based 
on  mean  free  path  compared  to  distance  from  the  source,  and 
we  are  looking  for  solutions  for  small  values  of  this  ratio. 
Once  we  have  the  fundamental  solution  for  small  Knudsen  number 
based  on  distance,  we  can  apply  it  to  problems  for  all 
Knudsen  numbers  based  on  mean  free  path  compared  to  body  size. 
Before  discussing  the  actual  nature  of  the  solutions,  we 
consider  the  application  of  the  fundamental  solution.   In 
general,  in  order  to  solve  for  flow  past  a  specific  body,  we 
cannot  avoid  the  travail  of  sxamming  sources  and  solving  the 
resulting  integral  equations.   If  we  are  successful  with  the 
mathematics,  a  detailed  picture  of  flow  is  given  outside  of  a 
region  of  several  mean  free  paths  surroiinding  the  body.   If  we 
choose  to  be  less  demanding  in  the  amount  of  information  that 
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we  want,  we  find  a  cache  of  information  directly  from  the 
fundamental  solution  itself.   It  is  shown  in  Section  5  that 
for  given  values  of  drag,  lift  and  heat  addition,  the  funda- 
mental solution,  with  these  values  inserted,  gives  the 
asymptotic  flow  field.  This  solution  is  asymptotic  in  the 
two  Knudsen  ratios  mentioned  in  the  preceeding  paragraph. 
For  objects  whose  dimensions  are  of  the  order  of  the  mean 
free  path  or  smaller  we  would  presume  that  the  fundamental 
solution  begins  to  describe  the  flow  many  mean  free  paths 
from  the  body.   For  objects  whose  dimensions  are  larger  than 
the  mean  free  path  (i.e.,  for  those  flows  which  embrace  fluid 
mechanics)  we  must  place  ourselves  many  body  lengths  from  the 
object  before  the  fundamental  solution  begins  to  describe  the 
flow  field.   In  either  case  if  the  problem  is  solved  by 
s\jmming  sources,  the  two  solutions  approach  each  other 
asymptotically  as  we  allcw  R  to  increase  (Sec.  3)-   We  should 
bear  in  mind  that  in  the  former  case,  that  of  using  the 
fundamental  solution  directly,  we  must  be  content  with  being 
ignorant  of  the  lift,  drag  and  heat  transfer,  whereas  in  the 
latter  case  the  lift,  drag  and  heat  transfer  is  computed  from 
the  calculated  flow  field. 

We  now  turn  to  the  approximate  fundamental  solutions 
which  are  determined  in  Sections  7  and  8.  The  analytical 
results  which  are  mentioned  refer  to  the  two-dimensional  case. 
This  is  the  case  dealt  with  in  this  report.   Let  us  note  that 
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since  the  terms  used  In  the  computation  come  from  the  Oseen 
grouping,  we  are  also  discuGsing  the  fundamental  solution 
of  the  Oseen  equations.  The  Mach  layer  solutions  of  Section 
8  are  uniformly  valid  with  respect  to  Mach  number,  as  we  pass 
from  supersonic  to  hi^h  subsonic  velocities.   Section  7  con- 
tains the  wake  calculation  which  is  uniformly  valid  for  all 
Mach  nxiJTibers.   Aside  from  the  wake  no  higher  approximation 
is  necessary  for  flows  up  to  high  subsonic  velocities.   For 
these  cases  the  inviscid  solutions  themselves  give  a  consistent 
picture  outside  the  wake.   In  discussing  the  results  of 
Sections  7  and  8,  we  deal  separately  with  transonic  and  super- 
sonic flow.   A  fluid  dynamical  approach  to  the  supersonic 
Oseen  equations  has  been  considered  by  Lagerstrom,  Cole, 
Trilling  [18],  also  Cole  and  Wu  [19].   Although  different 
methods  were  used,  there  is  agreement  in  the  results. 

For  supersonic  flow,  the  Mach  lines,  and  the  wake  line 
of  Inviscid  theory  act  as  a  skeleton  for  the  higher  order  flow 
field.  The  inclusion  of  higher  order  terms  fuszes  out  these 
lines.   We  find  that  the  diffusive  region  surrounding  the 

Mach  lines  is  centered  on  the  Mach  line  and  grows  with  the 

1/2 
iistance  along  the  Mach  line,  h,  as  (Lh)^   (Sf.c.  6  and  8). 

Furthermore,  the  fine  structure  of  this  layer  is  governed  by 

the  similarity  variable  (t'^/hL),  where  t  is  the  coordinate 

normal  to  the  Mach  line  (Sec.  6  and  8).   A  diffusion  operator 

is  responsible  for  the  Mach  layer,  and,  as  we  vjci.ld  expect, 

a  diffusion  operator  also  governs  the  wake  (Sec.  6).   The  wake 
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1/2 
spreads  as  (Lx)  '  ,  where  x  is  the  downstream  variable,  and 

the  fine  structure  of  the  v/ake  Is  governed  by  the  similarity 
variable  (y  /Lx),  y  being  the  transverse  variable  (Sec.  7). 
Allowing  the  mean  free  path  to  vanish,  we  recover  the  invlscld 
soluticns  (Sec-  4,  7,  8),   In  doing  this  we  go  from  a  para- 
bolic to  a  hyperbolic  partial  differential  equation  (Sec.  6), 
so  that  the  singular  lines  exist  as  unveiled  characteristics. 
A  most  important  feature  of  the  limiting  process  is  that  the 
convergence  is  not  uniform,  it  being  slower  more  distant  from 
the  origin. 

To  gain  Insight  into  this  non-uniform  convergence,  we 
look  at  it  in  another  but  equivalent  way.   In  looking  at  the 
supersonic  solution,  (Sec.  8)^,  and  also  from  remarks  above,  we 
see  that  h  and  L  are  interchangeable.  That  is,  allowing  h  to 
vanish  is  equivalent  to  a  vanishing  L.  But  as  we  mentioned 
before  (also  Sec,  8)  a  vanishing  L  gives  us  the  inviscia 
solutions,  and  also  therefore,  so  does  a  vanishing  h.   So  that 
for  a  small  but  fixed  L,  as  we  allow  h  to  approach  zero,  the 
solution  more  closely  resembles  the  invlscid  solutions.  This 
fortifies  the  idea  that  invlscld  flow  occurs  just  outside  the 
boundary  layer.   However,  no  matter  how  small  L  >  0  is  chosen, 
we  can  go  far  enough  away  from  the  origin  (and  hence  a  body), 
such  that  the  Mach  lines  are  diffused,  and  the  invlscld  solu- 
tion does  not  give  an  adequate  description  of  the  flow  there. 
We  are  now  in  a  position  to  give  an  over-all  picture  of 
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supersonic  flow  according  to  the  linearized  framework  we 
have  adopted.   Starting  from  the  body  we  first  have  a  boundary 
layer  region   (as  we  have  mentioned  earlier,  the  slip  flow 
region  cannot  be  described  by  our  analysis)  which  changes  to 
the  wake  as  we  go  doivnstream.   Outside  of  the  boundary  layer, 
the  flov;  becomes  inviscid  and  remains  so  for  distances  which 
are  0(3^  /L).   This  is  the  distance  at  which  the  layer  is 
approximately  double  the  body  size.   Sufficiently  far  away, 
the  region  in  which  disturbances  propagate,  i.e.,  the  Mach 
region,  and  the  wake  (which  grow  as  (Lh)  '     ,  (Lx)  ' 
respectively)  widen  enough  so  that  a  viscous  description  is   , 
again  needed.   With  the  exception  of  the  Mach  line  and  wake 
line,  any  ray  from  the  origin  ultimately  passes  into  and 
stays  in  a  region  of  inviscid  flow.   On  the  other  hand,  any 
parallel  to  Mach  lines  or  the  wake  line  ultimately  passes 
into  and  stays  in  a  viscous  region. 

A  similar  set  of  circumstances  also  holds  for  the  transonic 
case;  however  the  results  are  somewhat  more  striking.   Con- 
sideration of  the  wake  is  independent  of  Mach  number  and  the 
discussion  of  it  is  not  repeated.   For  simplicity  only  sonic 
flow  will  be  described.   Again  v;e  find  a  diffusive  layer 
centered  on  the  Mach  line,  which  in  this  case  is  the  y-axis. 
The  layer  thickness  grows  as  {-^  !•')'      (Sec.  8)  which  is  more 

Even  though  we  have  a  linearized  theory,  a  boundary  layer 
is  formed,  see  for  instance  [19]. 
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rapid  than  in  the  supersonic  case.   Detailed  information  of 
the  structure  of  the  layer  is  now  described  by  the  similarity 
variable  (xVy  L)  (Sec.  8).   The  differential  operator  which 
governs  the  layer  is  third  order  parabolic,  and  not  of  a 
well-known  form.   Perhaps  the  most  interesting  question  that 
can  be  asked  of  the  sonic  solution  is  its  limit  for  a 
vanishing  mean  free  path.   As  we  show  in  Section  8,  this  is 
the  much  maligned  solution  of  the  inviscid  sonic  equation 
(Sec.  4).  The  latter  violates  the  linearization  restriction, 
the  velocity  perburbation  grows  linearly  in  y,  and  predicts 
an  infinite  drag  because  of  it  (Sec.  4).   The  clue  to  the 
situation  lies  in  the  non-uniform  convergence  of  the  solution 
as  L  is  allowed  to  become  small.   No  matter  how  small  L  is 
made,  sufficiently  far  away  the  viscous  solution  takes  over 
and  decays  as  l/y  ^^ .      The  drag  then  is  bounded,  see  Section  8. 
In  the  last  two  sections  we  sketched  the  flows  for  super- 
sonic and  just  sonic  speeds.   Sections  7  and  8  give  the  uniform 
solution  from  which  the  sketches  were  taken.   It  is  exhibited 
in  closed  form  in  terms  of  an  Airy  function.  The  analysis 
of  Section  6  shows  that  the  solution  is  also  applicable  to 
high  subsonic  velocities. 
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3.   Governing  Equations. 

For  the  sake  of  preliminary  discussion,  we  shall  discuss 
the  steady  flow  about  an  arbitrary  body  by  means  of  kinetic 
theory.   As  a  governing  equation,  we  choose  the  single 
relaxation  time  kinetic  model  ([1],  [2]), 

(5.1)     4  •  Vf  =  v(fQ-  f), 

which  has  been  studied  by  Gross  and  his  co-workers  ([3],  [^]» 
[5]).   The  distribution  function,  f,  is  normalized  to  be  a 
mass  density.   The  collision  frequency,  v,  is  proportional 
to  the  mass  density,   f  denotes  the  locally  Maxwellian  dis- 
tribution function  described  below,  and  i   is  the  molecular 
velocity.  We  will  not  consider  the  approximations  needed 
to  arrive  at  (3-1),  since  this  is  discussed  in  the  literature 
{[J>] ,    [^])-   Besides,  before  very  long  we  shall  see  equation 

(3.1)  in  light  of  the  description  it  produces,  and  may  evalu- 
ate it  from  that  standpoint. 

H.  Grad  [6]  studied  the  fundamental  solution  of  (3-1)^ 
i.e.,  the  solution  of 

(3.2)  ^  •  Vf  =  v(fQ-  f)  +  (5(e)6(x). 

The  source  term  6{i) ,   will  be  discussed  presently.   Grad 
developed  a  technique  for  gaining  an  asymptotic  form  of  the 
fundamental  solution,  and  the  present  problem,  suggested  by 
him,  is  in  a  sense,  to  be  considered  the  next  step  in  the 
asymptotic  development. 
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The  single  relaxation  model,  simple  as  it  seems.  Is  still 
a  non-linear  integro-differential  equation.   This  is  so,  by- 
virtue  of  the  local  Maxwellian,  f  ,  which  contains  the  five 

o 

molecular  velocity  moments  of  the  distribution  function,  f , 
giving  the  local  density,  temperature,  and  velocity,  p(x), 
T(x)  and  u(x) . 

,[^-u(x)]2/2RT(x) 


(3.5) 

where 
(3.4) 


f  =  P(^)  ^. 

°   [27rRT(x)]^/^ 


1 
u 

RT 


1 


n 


v,_ 


As  is  shown  in  reference  [J>] ,   the  single  relaxation  model  can 
result  in  a  certain  way  from  the  weakly  linearized  Boltzmann 
equation;  however,  it  still  leads  to  a  non-linear  hydrodynamic 
theory,  as  can  be  seen  by  taking  velocity  moments  of  (5.1). 
Following  Grad  [6]  we  write 

""f  =  f°  +  g 


(3.5) 


p 

=  Po 

+ 

p 

T 

=    ^0 

+ 

T 

P 

=  Po 

+ 

P 

u 

+ 

u 

where 
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p^    -(4-u  )  /2RT 


is  the  unperturbed  constant  flow  at  infinity.   Using  the 
same  letter  for  a  quantity  and  its  perturbation  should  not 
lead  to  confusion  since  the  letter  will  henceforth  denote 
just  the  perturbation.   Denoting  the  relative  molecular 
velocity  by 


(3.7) 


c  =  i   -  u^ 


we  have 


(3.8)  { 


P  =  /  g  d^ 

pQ^  +  P^o  =  ^  ^S  d^ 

PqU  =  /  eg  d4 

p  =  i  /  c^S  d^ 

p^RT  =  /  (^  c^-  RTQ)g  d^ 

p  =  p  +  T 


Linearizing  the  local  Maxwellian  about  the  absolute  Maxwellian 

2 


(3.9) 


>  ^  f(0)  /i  .  £_  ,  /  c^    3x  T   ,  c-u  ) 

r   Po    ^2RT^    2^  T^   RT^  J  ^ 


and  introducing  the  above  into  equation  (5.1),  we  get 

(3.10)    i    .  7g  +  vg  =  vf(°^  l^-  +  (-^  -  3)  T_   c-u 

W    ;     ^   vfe  -t-  VB    vi     ]  p^  ^  ^2RT     2''  T   ^  RT 

■     (^^o       o       o     o 
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Next  Introducing  dlmenslonless  quantities 

r 


P    =    ^       ^  ^    = 


(3.11)     I 


u   = 


c    = 


o 
u 


'^   _  p 

T7    '      P  "  p; 


u 


^ 


,         u      = 


Vrt 


°   ,  X-^ 


VRT 


o 


g  =   (RTq)^/^  s/Pq 


1  2   ^ 

CO  =  ^7^  e 

^  (27r)^/2 


/RT" 


and  the  mean  free  path,  L,  by 


(3.11')   1=-^ 


o 


we  have  from  (3.8) 


P  =  J  S  d4 


(3.12)   . 


u  =  J  eg  d4 

P  =  /  i  c^g  df 


and  (3-1)  now  becomes 

As  a  result  of  the  linearization,  v  and  hence  L,  are  now 
constants.  This  is  the  equation  that  H.  Grad  ultimately 
studied  in  reference  [6]. 
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For  the  moment  we  make  the  spatial  variable,  x,  dimension- 
less  by  comparing  it  with  R,  which  is  representative  of  the 
distance  from  an  observer  point  to  the  origin. 

(3.14)  x/r  =  X'. 

With  this  stretching,  equation  (3.13)  becomes 

(3.15)  el  .  ^g  +  g  =  0)^^  +  (|c2-  |)T  +  S'.u 

where 

€  =  L/R. 

Taking  the  Fourier  transform  of  equation  (3-15)  and  denoting 
the  transformed  variable  by  the  untransformed  symbol,  for 
instance 

(3.16)  p(k)  =  /  p(x)  e"^^-^  d5? 
we  get 

(5.17)  gdc,?)  =  ^_p|j^[p  +  (1  c2-  |)T  +  c.s]  - 


e 

1  +  eik-^ 


^  ge-i^^'  ^ds 


J 


'2 

The  integral  is  over  the  surface  of  the  body,  2,  the  normal 

is  directed  into  the  body.   It  arises  by  use  of  the  divergence 
theorem,  from  the  parts  integration  of  the  Fourier  transform 


* 

The  surface  integral  can  be  taken  at  the  body  only  if  the 

object  is  slender,  and  the  boundary  conditions  allow  lineari- 
zation.  If  this  is  not  the  case,  the  surface  integral  is  to 
be  taken  at  some  surface  surrounding  the  body  at  which  the 
linearization  holds. 
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of  the  convectlve  term  of  the  differential  term.   The 
argument  of  the  exponent,  x',  Is  0(^/r)  where  ^  Is  a  length 
representative  of  the  length  of  the  body.   Expanding  the 
exponential,  we  get  an  asymptotic  representation  of  the 
Integral  in  the  small  parameter  (^/r) 

(5.18)    /  ge"^^^'  |.ds  -^  /  g^ds  -  /  Ik-x'  g4-ds  +  ... 
2  2         2 

The  leading  term  of  equation  (3-18)  Is  0(<^^/r^),  by  virtue 

of  being  a  surface  Integral.   At  this  point  we  have  two  small 

dlmensionless  ratios,  l/r  and=2!^/R.   There  Is  no  reason  to 

compare  them  in  the  analysis.  The  first  ratio  e  =  l/r,  as 

we  will  see  shortly  after  expanding  in  powers  of  it,  acts  as 

a  guide  in  choosing  approximate  equations;  whereas  the  second 

ratio  v^f/R  as  we  see  from  equation  (3.I8)  splits  the  boundary 

conditions  into  orders  of  ^/r.   In  what  follows  we  shall 

carry  only  the  first  term  of  the  expansion  (3.18).   This 

then  yields  a  solution  which  is  asymptotic  in  the  ratio  .^/r. 

By  virtue  of  the  forthcoming  expansion  in  e,  this  solution 

is  asymptotic  in  the  two  ratios  ^f/R  and  l/r.   The  expansion 

in  €  =  l/r  demands  that  the  solution  be  used  only  many  mean 

free  paths  away. 

Returning  to  equation  (3.I7),  retaining  only  the  first 

term  of  the  integral,  (5.18),  and  removing  the  stretching 

(3-1^),  we  get 

D.19)       g(k,4)  =  — ^^-^  (p  +(i  c2-  |)T  +  c-uU  -Mli4 

l+lLk-4     L       -  ^  -^       l+lLk-^ 
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with 

(3.20)    6{i)   =  -  /  g(x,fK-<is. 

2 

Henceforth,  when  expanding  with  respect  to  the  mean  free 
path  L,  we  shall  mean  that  L  is  small  in  comparison  to  the 
radial  distance  of  an  observer  point.   From  (5.20)  we  see 
that  6{i)   is  the  mass  flow  out  of  the  body  per  volume  of 
velocity  space  (this  is  not  dimensionally  true,  due  to  the 
dimensionalization  (3.11)).   In  addition  we  see  that  (3. 19) 
is  precisely  the  equation  we  would  get  by  our  procedure  if 
we  had  considered  Grad's  equation,  (3.2),  to  begin  with. 
We  therefore,  see  that  solution  of  (3-19)  applies  to  two 
different  cases.   First  it  gives  the  asymptotic  solution  for 
flow  past  a  finite  object  subject  to  the  preceding  remarks. 
Second,  it  gives  a  fundamental  solution  for  the  strongly 
linearized  single  relaxation  kinetic  model.  This  distinction 
is  not  just  a  matter  of  interpretation.   Since  we  have  linearized 
our  equation,  the  solution,  when  considered  as  applied  to  the 
asymptotic  gas  flow  past  an  object,  must  agree  with  the 
linearization.   Hence  it  applies  only  to  those  regions  where 
products  of  the  hydrodynamical  quantities  are  small  compared 
to  the  quantities  themselves.  This  restriction  is  not  to  be 
invoked  on  the  solution  when  it  is  viewed  as  a  fundamental 
solution.  The  fundamental  solution  is  used  to  construct  solu- 
tions of  boundary  value  problems,  and  it  is  only  the  latter 
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solutions  which  must  comply  with  the  linearization  requirement, 
Furthermore,  the  latter  case  places  no  restrictions  on  the 
dimensions  of  the  object.   In  fact  we  can  use  the  fundamental 
solution  to  construct  solutions  of  boundary  value  problems 
through  all  ranges  of  Knudsen  number. 

On  multiplying  equation  (3.19)  by  1,  c,  (i- 'c  -  1)  and 
integrating  in  velocity  space  we  get  an  inhomogeneous  linear 
system  of  equations  for  p,  T,   u,  which  Grad  [6]  first  derived. 
The  equations  are 


(5.21) 


with 


1  -  /  Qdt 


-I   I  (c2-5)fid4    -/  c^nde 


-/  i(^'2-3)i^dt  1-/  ^{^^-3fm     -j\  ^^{c'^-J>)m 


■J   c^i7d4 


-/i-  c^(c2-3)fid?  1-!  Z^c^m 


2  r' 


/^i 


=  L 


6r 


V 


(3.22)  I 


p  - 

0) 

l+iLk-4 

6       -    f        ^^^ 

^       J    1+iLk-^ 

1 

6       - 

r  (1/3  c^-Dsdl' 

2        J 

1+iLk-? 

S'\ 

r     c  d  d| 

1  +i  T.k- .  ^ 

L. 
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The  matrix  equation  (3. 21)  is  a  5X5  system,  repeated  sub- 
scripts representing  summation,  and  unrepeated  subscripts 
ranging  through  the  number  of  dimensions  of  configuration 
space. 

In  order  to  perform  the  indicated  integrations,  we 
expand 

(3-23)    fi  ^CD  [1  -  lUc'i   +  L^{k''()^   +  ,  .. 

The  resulting  Integrated  series  gives  the  asymptotic 
representation  of  the  integrals  occuring  in  (3.21).   Retaining 
only  the  first  term,  i.e.,  dropping  0(L),  equation  (3.21) 
Identically  vanishes.   Carrying  terms  of  0(L),  i.e.,  taking 

(3-24)    fi  --  (D  (1  -  ILk-f) 


we  get 


(3^25) 


/iu^.k    0 


V 


Ik. 


lu^'k 
o 


Ik. 


iu^-k6^^ 
o   rs 


This  is  the  system  studied  by  H.  Grad  in  his  report  [6].   It 
is  simple  to  see  by  direct  Fourier  inversion,  that  the  system 
(3.25),  is  just  the  linearized  Euler  equations  with  source 
terms.   The  results  that  Grad  derived  are  essentially  those 
usually  gotten  from  Inviscid  flow  theory.   As  is  well-known 
from  gas  dynamics,  the  solution  of  equation  (3.25)  is 
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singular  along  the  Mach  cone  for  supersonic  flow  (see  Sec.  4). 
H.  Grad  [6],  also  showed  the  existence  of  a  singularity  along 
the  flow  axis  behind  the  origin  independently  of  Mach  number. 
The  divergence  arises  from  the  fact  that  Grad  solved  equations 
(3.25),  whereas  the  comparable  potential  equation  is  solved 
in  gas  dynamics.   The  additional  singularity  arises  by  not 
excluding  rotational  flow  (as  is  implicit  in  the  potential 
equation),  and  is  the  inviscld  of  limit  of  the  wake  behind 
the  body.   Grad's  solution  is  the  correct  fundamental  solution 
of  the  Euler  equations  without  a  mass  source. 


4.   Solutions  to  the  Euler  Approximation. 

Before  going  on  to  the  higher  approximations,  we  work 
out  the  two  dimensional  analogue  to  Grad's  analysis. 

In  order  to  carry  out  the  inversion  necessary  for  the 
solution  of  (5-25),  we  must  first  consider  the  source  terms 
which  occur  on  the  right  hand  side  of  the  equation.   As  was 
shown  by  H.  Grad  in  [6],  the  Fourier  inversion  of  a  term 
such  as  G,  appears  as  an  object  of  size  0(L)  located  at  the 
origin.   Physically  this  means  that  the  source  appears 
smeared  out  to  the  extent  of  a  mean  free  path.   This  smeared 
source  can  also  be  expanded  as  was  the  integral  term  in 
equation  (5.I7).   Actually,  this  operation  can  be  carried  out 
earlier  by  expanding  6  in  the  same  way  as  indicated  in  (3.24) 
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This  leads  to  an  asymptotic  representation.  In  the  sense  of 
distributions,  and  is  just  a  multipole  expansion  in  terms  of 
the  Dirac  5-function  and  its  derivatives.  Therefore,  we  have 
to  the  lowest  order 

f 


(^.1) 


6^  ^  /  dd^ 


V- 


6-,^    J   c6dC 


For  the  sake  of  simplicity  we  take 


(4.2) 


/   6(Od4   =   0 


which  is  a  statement  of  mass  conservation,  and 


(^.3) 


/-< .  .-^2  "a 


/  6(orcae  =  0 


which  states  that  no  net  energy  exchange  is  taking  place 
between  the  body  and  the  fluid.   We  shall  later  show,  in 
Section  9,   the  effect  of  heat  addition.   With  these  condi- 
tions, we  take 


r 


(4.4) 


d-,  ^  TU^ 

5     o 


2  ,  2 
5   o 


Therefore,  t  is  a  negative  constant  representative  of  the 
drag,  and  is  of  the  dimensions  of  length.  These  conservation 
statements  are  only  sensible  in  talking  about  the  asymptotic 
flow  past  a  body.   In  talking  about  the  fundamental  solution 
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we  can  still  make  the  statement  about  the  conservation  of 
mass,  as  long  as  we  agree  not  to  consider  flows  with  mass 
addition.   On  the  other  hand  we  no  longer  have  the  right 
to  avail  ourselves  of  the  energy  equillbriiim  equation  (4.5), 
and  the  exclusion  of  lift,  as  is  implicit  in  equation  (4.4). 
These  are  in  general  needed  when  one  attempts  to  actually 
construct  solutions.   We  will  return  to  these  remarks  in 
Section  9*   For  the  present  we  assume  equations  (4.5)  and 
(4.4)  for  the  sake  of  simplicity. 

We  can  simplify  our  calculation  by  noticing  the  4X4 
matrix  (5*25)  can  be  reduced  to  a  5X5  matrix,  if  we  solve 
for  p,  T,  and  k'U.   This  is  accomplished  by  taking  the  inner 
product  of  the  last  row  of  equation  (5.25)  with  k.   In  this 
case  (5«25)  becomes 


(4.5) 


i(u^-k) 


ik' 


i(u^-k) 


ik' 


1      \ 
i(u^.k) 


/ 


T 


k«u 


2-^2 
5   o 

\  Tk'ul 


and  inverting  the  matrix 


(4.6) 


P 


T 


/|k2-(a^.k)2     -k^ 


k-u 


1 
D 


-  2^2 


(u^-k)k^ 


u  -k 
o 


k'^-(u^-k)''  4  u^'k 
^  o   '    50 

(u^.k)k2   -(u^-k) 


with 
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(4,7)  D  =  i(u^.k)[|k2-   (S'o-k)2] 


=   ^^0^1    t(f  -u,2)k,2^f  k/] 


3 

i\fhere  the  upstream  velocity,   u    .   is  taken  as  oriented  along 
the  X, -axis.      Since 

(4.8)  u  ^  =  I  M^ 

Op 

2 
where  M  is  the  Mach  number,  the  coefficient  of  k,   in  D  is 

o 
proportional  to  (1-M  ).  Because  the  character  of  the  equations 

changes  from  elliptic  to  hyperbolic  as  M  increases  through 

the  value  one,  we  consider  each  regime  separately. 


Supersonic  Case 

(4.9)     (1  -  |u/)  =--^  ,  P^  >  0. 

Solving  for  the  density  from  the  matrix  equation  (4.6) 

2  ,5  „  2  ^  2  ,  2v  1 


(4.10)    ^  =  TU^^  (^k^-^^  +f  V)  fe 


and  taking  a  Fourier  inversion 

-3i  TiT  p^  ^   v2 
(4.11)    p  =  _-2!_  (IL 


■^^x2 


o         #  ik, x+ik^y 
+  2  Sl_)   1   I  e  -^    ^  dk 


3  ^„2'  /^_s2 


5y   (2Tr)^J  k^(k^-pk2)(k3^+pk2) 
Similarly  we  get  for  the  temperature  T 
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(4.12) 


T  = 


2i  TU  p' 


[-  u 


2  S' 


^4. 


°  3x^ 


ik,x+lk:oy  

e  -^    ^  dk 


Sy^   (27r)^   k^(k^-pk2)(k^+pk2) 

Finally,  in  order  to  determine  the  velocity  field,  we 
return  to  (5«25)  to  get 


(^.15)    u^  = 


-k^(p  +  T)   ItUq 


%^1 


Kl 


with  u  ,  a  unit  vector.   Inversion  gives 


+iT(5P  +2)     :. 
(4.14)    u  =  1^ -V     ^ 


r  ^ 


(27r)' 


r 


J 


ik-iX+ik^y  

e  -^    ^  dk 

k-^(k^-3k2)(k^+pk2) 


u. 


-  Ti 


(27r)' 


ik,x+ik^y 
e   -^    ^  dk 
kn 


In  the  evaluation  of  the  two  integrals  which  occur,  the  path 
of  integration  in  the  complex  k-|-plane  is  taken  slightly 
below  the  reaT  .xis.  This  eliminates  advanced  solutions, 
i.e.,  solutions  propagating  upstream.  The  integrations  are 
easily  carried  out,  in  the  sense  of  distributions,  see 
Lighthill  [17],  and  we  get 


r 


(4.15) 


(27r)' 


ik'X 
e     3x 


=  iH(x)5(y) 


J 


where  H(x)  is  the  Heaviside  step-function,  and  6(y)  the  Dirac 
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delta-function.  The  second  Integral  Is 

r 

1 
(27r)^   k. 


(^.16) 


J 


e H =iH(|l(ly| 


-  -  |y+Px 


On  substitution  Into  equations  (4.11),  (4.12)  and  (4.l4) 
we  get 

(4.17)    P  =  f  TU^6(y)H(x)  -  lo'u^cL  ^  (^u^^  5)[6(y  cos  0  + 

-^  o 

+  X  sin  0)  +  5(y  cos  0  -  x  sin  0)] 


rJ 


(4.18)    T  =  -  I  T{^^6(y)H(x)-.^3^g(g^^g  (2u2+  5)[6(y  cos  0  + 


+  X  sin  0)  +  5(y  cos  0  -  x  sin  0)] 


(^.19) 


Un 


u. 


0 


TH(x)6(y)  + 


tH(x)(2u>  5)  /sec  0 


6u, 


CSC  0 


5(y  cos  0  +  X  Sin  0)  + 


tH(x)(2u^+  5)  /sec  0  \  5( 
6u^       \^-csc  0  y 


y  cos  9  -  X  sin  0) 


X 


Fig.  4.1 
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The  solution  exhibited  in  equations  (4.17-4.18)  are 
non-zero  only  on  the  heavily  marked  lines  in  Figure  4.1. 
With  the  exception  of  the  line  concentration  5(y)H(x),  the 
solutions  are  what  we  would  anticipate  from  gas  dynamic 
potential  theory.   A  discussion  of  these  singularities  in 
the  three  dimensional  can  be  found  in  [6],  (the  only 
qualitative  difference  in  the  solutions  being  that  the 
three  dimensional  supersonic  case  exhibits  debris  behind 
the  Mach  cone,  corresponding  to  the  tail  of  the  solution  to 
the  wave  equation  in  two  space  dimensions).   It  suffices  here 
to  say  that,  when  the  supersonic  flow  past  an  object  is 
viewed  in  the  inviscid  limit,  the  viscous  layer,  i.e.,  the 
Mach  lines  and  the  wake,  are  compressed  into  line  discontinu- 
ities.  It  is  one  of  the  goals  of  this  report  to  examine  these 
regions. 

Subsonic  Flow:   (1  -  ^  u^  )  =  i^  ,  a^  >  0. 

The  subsonic  case  is  found  in  essentially  the  same  way 
as  the  supersonic  case.   Rather  than  repeat  the  analysis,  we 
just  give  the  salient  forms  which  appear. 

After  solving  for  the  transformed  flow  field  and 
inverting  the  Fourier  transforms,  we  get 


P  = 


^  _     o  /5  S    .  2  S"  N   a    I  e  "    "  dk 
—  w  ■  o  ^   ^  o'' 


5    5  Sx^   5  By^  (27r)^  1  k^(k^-iak2)(k3^+iak2) 
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(4.20) 


T  =  — - —  (u 


r 


r)-^ 


°   ax^        Sy2    (2Tr)2 


e     -^  ^     dk 

k,  (k, -lakp )  (k,  H-iakg ) 


u 


1\  =  ^    (5a"-2) 


u. 


lk,x+ik^y  

e     -^  "^     dk 


-  It 


S     j  {2it)'^    I    (k^-iak2)(k-L+iak2) 


ik^x+lk^y  

e      -^  "^     dk 


0/   (27r)' 


As  before  the  Integrals,   in  the  ordinary  sense,   are  divergent, 
"but  exist  In  the   sense  of  distributions,    [17]  • 


(^.21) 


f    lk-,x  +  Ik^y 
e     -^  "^     dk 


(27r)' 


=  15(y)H{x) 


(4.22) 


^    lk,x  +  Ik^y  

1  e     -^  ^     dk 


(2Tr)^J    (k^-lak2)(k3^+lak2)  ^^^ 


log  r 


(4.25) 


2  r       i^l^     +     '-^2^    -rr- 

a  !    e  dk 


{2Tr  J   k^(k^-lak2)(k^+lak2) 


^   _  i|y|H(x)    _ 


with 


lo^logr+iy  tan-^   (^) 
27r  2Tr  ^ox' 


(4.24) 


=    Va^x^  +  y^ 


Inserting  these  into  the  flow  field  (4.20),  we  have 
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r 


(^.25)  { 


P 


T 


2tu 


TU. 


-H(x)6(y)  +10^^  [5a   -   2] 


Trr 


2tu 


^5(y)H(x)  +1-^   [5a2-  2] 


Trr' 


log  r 


The  subsonic  case  given  above  is  listed  only  for  completeness. 
With  the  exception  of  the  wake  and  those  high  subsonic 
velocities  which  may  be  considered  transonic,  it  will  not  be 
studied. 


Sonic  Flov/:  u 


o 


5 
3 


Although  the  linearized  inviscid  equations  give  seemingly 
absurd  results  for  just  sonic  flow,  it  is  of  interest  for 
comparison  with  later  work  to  have  this  computation.   Returning 
to  equations  (4.11),  (4.12)  and  (4.14)  we  have  in  this  case: 


(4.26) 


(^.27) 


(4.28) 


_  ___  X  2  5"^  N   1     e  -^      ^ 


p  =  ^  (5  ^ 


dk 


^l'^2 


_  2  iT  ,5  ^ 


T  =  ^ 


(?^- 


■) 


r  ik^x  +  ik^y  


"r  = 


^  %  ^  ax2   ^y2'  \2^ 

f    ik,x  +  ik^ 
_  t5  r^    1    1  e  ^      ^ 


dk 


k^k^^ 


y 


V, 


dk 


5  "(i^O^ 


iTU^ 
O 

(27r)2j 


/-  Ikx  +  iky 


dk 


kn 
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The  Integrations  which  occur  are  straightforward,  as  before 


(^.29) 


ik,x  +  lk„y  

1    I  e  ^      l_^  ^   _.jj(^j  jyj 


(27r)' 


k^kg 


r  Ik^x  +  ik^y 


(4.30) 


(Stt)' 


dk 


=  -  5(x) 


J 


Substituting  into  equations  (4.22-4.24)  we  get 


r 


(4.51) 


P  =  2^  (|  5.(x)|y|  +|H(x)5(y)) 


T  =-^  (4  6'(x)|y|  -  |H(x)6(y)) 
5Uo  5  ^ 


u. 


u. 


v_ 


+1    |TH(x)5(y) 


This  solution  grows  linearily  in  the  y-direction. 
Because  of  this,  the  assumption  of  linearization  soon  fails 
and  in  addition  it  is  simple  to  show  that  an  infinite  drag 
results.   Nevertheless,  these  solutions  become  meaningful 
when  we  examine  their  viscous  counterpart  in  Section  8. 
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The  Higher  Approximations. 


In  this  section  we  consider  the  system  of  equations  for 
u,  when  we  ta] 
by  equation  (3.22), 


p,  T,  u,  when  we  take  the  next  approximation  to  Cl,   defined 


(5.1)     (2  =  o)  'j  1  -  ILk-l  +  (iLk-O^J 


=  CD  [1 


ILk-c  -  IL(k.ir^)-  L^(k-c')^-  L^(k.S')k.UQ 


-  L' 


■c^-^o)'] 


When  this  Is  substituted  Into  the  matrix  equation  (3.21)  we 
get  a  linear,  Inhomogeneous  matrix  system,  as  before.  In  p, 
T,  u. 


(5.2) 


lu^.k+L(u^.k)' 
+  Lk^ 


2         2 
3 


Lk' 


ik_+2Lu_-kk. 


i(u^.k)+L(u^.k)2     I  (ik^  + 


+   7Lk2 


2Lu'  •  kk^ ) 
o        s 


lk^+2L(u^.k)k^        lk^+2L(u^.k)k^      [ia^-k+(u^.k)''L 

+k-L]5^3+2Lk^k3/     \ 


The  Integrations  Involved  in  determining  the  matrix  elements 
are  well-known  and  a  complete  list  Is  given  In  [12].   In 
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this  form  the  equations  of  (5' 2)  do  not  seem  to  bear  any 
relation  to  known  hydrodynamic  systems;  however,  this  is 
just  mldleadlng.   The  system  is  the  Navier-Stokes  set 
linearized  away  from  the  constant  free  stream  velocity. 
This  is  usually  referred  to  as  the  Oseen  equations.   For 
instance,  on  multiplying  the  system  (5-2)  by 


(50) 


1  +  iL(u  -k)    0    iLk„ 
^  o  r 


0 


we  get 


i(u^-k) 


(5.M 


2    2 


\ 


0 


ik. 


i(u^.k)+L(u  -kf  %   (i+2Lu  .k)k, 

^  O        ^  O  3  ^        OS 

+  lLk2 

l+2L(u^.k)]k^   (l+2L(uQ-k)]k^  [iu^-k+(u^-k)2L+ 


+1^  Ll5^^+2Lk  k 


rs    r  s 


u. 


o(l2)\ 


6, 


'3    / 


The  source  relations  of  Section  5  have  been  applied  to  the 
first  row.   From  the  first  row,  we  see  that  linearized 
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0 

iLk^                  \ 

l+L(u^-k)l 

f  ILk^ 

iLk^ 

l+lL(u^.k) 

o 
continuity  equation  is  satisfied  to  0(L  ).   The  linearized 

continuity  equation  is  in  fact  satisfied  identically,  as 

can  be  seen  by  merely  integrating  equation  {J>.2),   and  then 

linearizing.   To  get  the  familiar  Navier-Stokes  system  we 

have  only  to  multiply  (5.2)  by  the  matrix 


1  +  iL(u^-k) 
(5.5)   /  0 

iLk„ 


However,  since  there  is  no  compelling  reason  to  use  the  Navier- 
Stokes  system,  we  will  take  (5'^)  as  our  system. 

That  we  arrive  at  the  Navier-Stokes  set  should  cause 
us  no  great  surprise.   Without  going  into  details,  it  is 
simple  to  see  that  the  expansion  of  Cl   in  the  above  analysis 
is  identical  to  an  Iterative  prodecure.   This  mimics  the 
Chapman-Enskog  method  in  the  first  two  steps  (the  first  two 
steps  in  the  Chapman-Enskog  method  lead  to  the  Navier-Stokes 
set,  see  [16]). 

Since  we  have  a  system  equivalent  to  the  Navier-Stokes 
set,  we  have  also  a  calculation  of  the  transport  coefficients. 
These  will,  of  course,  depend  upon  the  choice  of  v.   A 
systematic  discussion  of  the  Boltzmann  equation  leading  to. 
the  single  relaxation  model,  is  given  by  Gross  and  Jackson   v 
[4].   The  value  of  v  that  the  analysis  indicates,  results  in 
a  viscosity  coefficient  which  Chapman  and  Cowling  call  the 
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first  approximation  to  the  viscosity  coefficient.  However, 
the  heat  conductivity'-  is  off  by  a  factor  of  3/?.  This  fact 
has  already  been  noticed  by  Gross  [5]. 


6.  Boundary  Layers. 

Starting  with  the  matrix  equation  (5-^)^  one  can  in 
principle  solve  for  the  flow  field  (p,  T,  u),  (indeed  even 
starting  with  (3.21)  the  method  of  solution  is  clear).   One 
has  only  to  invert  the  matrix  equation,  to  solve  for  the 
Fourier  transformed  field,  and  then  invert  the  transforms. 
Formally 


'-'  \ 


(6.1) 


\ 


p  \ 

T      ,  —    - 

1 

^ik-x 
e 

C 

■L       1 

ro_\n 

D 

u  / 


0 
6. 


dk 


where  n  is  the  dimension  of  the  space,  D  the  determinant  of 
the  matrix,  and  C  the  transpose  of  the  cofactor  matrix. 
However,  such  a  plan  does  not  seem  feasible.   Instead  we 
adopt  the  philosophy,  that  the   invlscid  solution  developed 
in  Section  3  is  adequate  except  in  the  neighborhood  of  the 
singular  lines.   That  is,  when  the  solution  of  the  invlscid 
equations  has  meaning  as  a  function  in  the  ordinary  sense 
it  is  adequate,  and  v;hen  the  solution  has  meaning  only  in  the 
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sense  of  distributions  It  points  out  an  inadequacy  In  the 
solution,  and  hence  in  the  governing  equation  in  that  region 
of  space.   In  the  remainder  of  this  section,  we  consider  the 
form  of  the  second  approximation,  Navier-Stokes,  in  the 
regions  of  singularity  of  the  first  approximation,  Euler. 

As  in  the  Euler  approximation  we  can  separate  our  system 
so  that  we  only  consider  p,  T,  and  k'U.  Taking  the  divergence 
of  the  last  row  of  i5-^)  >   the  matrix  equation  becomes 


(6.2) 


i(G^-k) 


|Lk2 


i(u^-k)+L(u^.k)2       I  (i+2Lu^.k) 


5 


[i+2LuQ-k]k^      {l+2L(uQ-k)^k' 


lu^.k+(uQ.k)^L 


+  3k^L 


P 


k-u 


/\  / 


^2 
k.6 


'/ 


The  determinant  of  the  matrix  is  easily  computed,  and  to 
0(L)  is  given  by 


(6.3) 


D  =  -i(u  .k)^+  1  I  (u^-k)k2+  %   k^L  -  2(u^-k)^L  + 


3  ^  o 


•  \    1^        \2  2 
+  ^  (uQ-k)VL. 
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As  can  be  seen  from  (6.1),  D,  if  Fourier  Inverted,  i.e., 
(6.4)     D  =  (i5„-V)^-  I  (Q^-VJV^  +  f  LV'-  2L(S^-V)'' 

+  I  (C„-V)V, 

is  the  differential  operator  governing  each  of  the  dependent 
variables  p,  T,  V-u.   Therefore,  we  can  expect  to  find  the 
correct  approximations  in  the  neighbpthood  of  the  singularities 
of  the  Euler  system,  by  examining  the  operator  D. 

Henceforth  we  restrict  attention  to  the  two-dimensional 
case,  and  we  will  Introduce  pertinent  results  of  the  three- 
dimensional  case  only  for  the  sake  of  comparison.   Let  us 
first  consider  the  neighborhood  of  the  Mach  line  as  in 
Figure  4.1.   First  we  rotate  our  coordinate  frame  through 
the  Mach  angle  Q,   as  in  Figure  6.1,  so  that  h  denotes  distance 
along  the  Mach  line,  and  t  the  coordinate  normal  to  the  Mach 


Fig.  6.1 


line.   We  assiome  that  the  Mach  line  singularity  is  actually  a 
narrow  layer  through  which  the  physical  variables  change,  and 
we  spread  the  layer  in  the  lateral  direction  by  factor  l/5, 
where  5  is  the  layer  thickness. 
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(6.5)  t/6  =  t' . 

This  should  actually  be  thought  of  as  -ry^^  5  =  VR,  where  R 
Is  the  large  distance  from  the  object  mentioned  in  Section 
3,  and  A  is  the  actual  thickness  of  the  layer.   Equivalent 
to  the  spreading  is  the  statement  that 

(6.6)  SM  -- 0(1/6), 

or  if  we  deal  with  wave  number  space,  that  the  wave  number 
corresponding  to  t  is  0(l/6).  We  now  carry  out  the  analysis 
in  wave  number  space,  i.e.,  we  consider  the  determinant  (6.3). 
Denoting  the  wave  niimber  vector  k  by 

(6.7)  k  =  (k^,  kg) 

we  perform  the  rotation  indicated  in  Figure  6.1,  and  represent 
it  by 


(6.8) 


/  k-,  \    /  -sin  6   cos  0\  /  m 
y  kp  /    \   cos  6   sin  9  /  \  n 


so  that  m  denotes  the  wave  number  across  the  layer,  and  n 
the  wave  number  along  the  layer.   Substituting  the  transforma- 
tion (6.8)  in  the  determinant  (6.5)  we  have 

(6.9)     D  =  -  I  iu  mn^sin  0  -  ^  lu  m^n  cos  e  -  ^  Lm^  + 

+  0(Lm^n,  Lm^n^,  Lmn',  Ln  ,  n^). 
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Use  has  been  made  of  the  relation 

(6.10)  u^^  sln^e  =  5/:5. 

Details  of  the  remainder  terms  In  (6.9)  are  not  given  since 
they  are  negligible  compared  to  the  terms  carried. 

Let  us  assume  that  9   is  neither  near  zero  nor  ir/2;   then 
cos  0,    sin  0'-  0{l).   Then  according  to  (6.6) 

(6.11)  -  -I  Lu^mn^sln  0   -^0(l/6) 


10 


(6.12)     -  ^  Lu^m^n  cos  6  ^  0(l/5'^) 


(6.13)     -  ^  Lm^-^  0(L/5^). 


Since  5  is  considered  to  be  a  small  parameter,  these  relations 

indicate  that  (6.12)  can  be  neglected  relative  to  (6.11). 

Furthermore,  since  the  Euler  term  is  inadequate  for  describing 

the  flow  in  the  region  of  Interest,  we  take  (6.11)  and  (6.I5) 

to  be  the  same  order  of  magnitude.   We  conclude  therefore 

that 

(6.14)    D-  -  ^  iu  m^n  cos  e  -  ¥  Lm^. 

Also  since  0(l/6  )  =  0(l/5  )  we  can  further  conclude  that  the 

1/2 
thickness  of  the  layer.  A,  varies  as  -^^  (RL)    .   Furthermore, 

from  (6.14)  we  see  that  the  layer  is  governed  by  a  diffusion 

phenomenon. 
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As  a  partial  justification  of  the  procedure  (fuller 
justification  comes  when  we  find  the  solution  is  in  agree- 
ment with  the  assumptions  made),  we  should  verify  that  the 

2 
next  order  terms,  i.e.,  L  terms,  are  Ignorable.   Terms  of 

o 
0(L  )  arise  in  two  different  ways,  first  those  terms  which  we 

have  been  neglecting  in  expanding  the  determinant  of  the 

matrix  in  equation  (6.2).   By  inspection  these  are  0(L  /^    ). 

Second,  those  which  arise  in  carrying  the  expansion  at  Q. 

one  step  further,  that  is,  taking 

(6.15)    ^^   1  -  (iLk'O  +  (iLk-U^  -  (iLk.^)^. 

We  shall  refer  to  the  system  of  equations  implied  by  (6.I5) 
as  "Burnett",  since  they  are  in  a  formal  sense  similar  to 
the  Burnett  equations.   Considering  this  system  we  again  find 
that  the  set  can  be  separated,  and  further  that  the  largest 
order  contribution  to  the  determinant  is  0(L  /5^). 

We  now  have  three  types  of  terms,  which  we  refer  to  as 
Euler,  Navler-Stokes  and  Biirnett;  and  which  we  estimate  to 
be  0(1/5^),  0(l/5  ),  0(L^/55),  respectively.   Assuming  that 
the  Euler  and  Navler-Stokes  terms  are  of  the  same  order, 
leads  to  a  6^0(L  '^  ).   If  this  is  substituted  into  the 
Burnett  term,  0(L  /^    ) ,    it  is  clear  that  the  Burnett  term  is 
negligible  under  this  assumption.   The  assumption  that  the 
Euler  and  Navler-Stokes  term.s  are  of  the  same  order  will  in  a 
sense  be  removed.   It  will  now  be  shown  to  be  the  only 
possibility. 
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As  we  have  mentioned  before,  we  expect  the  solution  to 
the  Euler  equations  to  give  an  adequate  description  of  the 
flow  In  the  greater  part  of  the  field.   Only  In  certain  narrow 
regions  do  the  Euler  solutions  break  down.   As  we  enter  these 
"narrow"  regions  from  the  Euler  region,  we  expect  on  physical 
grounds  (and  assume  mathematically)  that  the  differential 
operator  governing  this  region  changes  continuously.   In  other 
words,  there  are  differential  terms,  previously  neglected, 
which  gradually  become  as  Important  as  the  Euler  terms.   There 
Is  no  sudden  change  of  one  operator  into  another.   The  changing 
of  the  Euler  operator  into  the  Navier-Stokes  operator  has 
been  shown  in  the  last  paragraph  to  be  consistent.   If  we 
attempt  to  change  the  Euler  operator,  0(l/5  ),  into  the 
Biirnett  operator,  0{L^/b^) ,   we  find  that  6—  0(L^'^^).   This 
is  contradictory  since  it  implies  that  the  Navier-Stokes 
operator  is  of  a  lower  order. 

We  should  not  discount  the  possibility  that  after  the 
Euler  operator  has  changed  into  the  Navier-Stokes  operator, 
that  the  Navier-Stokes  operator  will  change  to  Burnett.   Equat- 
ing the  orders  of  the  Burnett  and  Navier-Stokes  term,  we  find 
that  such  a  transition  can  indeed  take  place  (the  Euler  terms 
are  of  higher  order),  but  in  a  sublayer  of  actual  thickness 
5  =  0(L).   It  seems  likely  however  that  in  such  a  region  the 
remainder  terms,  i.e.,  the  higher  order  terms  of  the  expansion 
in  L,  would  be  equally  Important.   In  other  words,  the  detail 
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which  the  Burnett  equations  try  to  furnish  In  the  sublayer 
can  only  be  gotten  from  the  full  kinetic  equation  itself. 

So  far  we  have  discussed  in  detail  the  approximations 
involved  in  the  purely  supersonic  case;  the  same  arguments 
apply  in  the  other  ranges.   In  the  almost  sonic  region, 
I  cos  9 1  «  0(5),  when  the  upstream  velocity  varies  from 
high  subsonic  to  low  supersonic  speeds,  the  operator  (6.9) 
is  approximately  given  by 

5  .^   2   10  T  2 


(6.16)    D  ^  -  ^  iu^mn""  -  ^  Lm^. 


All  arguments  presented  above  for  the  purely  supersonic 
case  apply  again.   Here,  however,  the  thickness  of  the  layer 
varies  as  ^  {'^'L)    '^ .      In  the  event  that  cos  0-->O(6),  it 
is  clear  that  all  three  terms  occuring  in  (6.9)  are  of  equal 
importance.   Therefore,  in  order  to  consider  a  solution  which 
is  uniformly  accurate  with  respect  to  Mach  number  in  all  the 
regions,  we  must  consider  the  entire  operator  D,  less  the 
remainder  term,  i.e.,  we  take 

(6.17)  D  =  -  5-  iu  mn  sin  0  -  -5-  iu  m  n  cos  Q   -   ^  Lm  . 
The  corresponding  differential  operator  is 

(6.18)  D  =  11^  ({T^sin  e  ^  ^r   2{^^cos  ^  ^  -  2L  ^) . 
In  the  supersonic  case  the  differential  operator  becomes 

(6.19)  x>.^^t^(^^cone\^-^t^) 
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and  as  indicated  before  the  layer  spreading  is  governed  by  a 
diffusion  equation.   In  this  case  we  expect  the  solution  to 
depend  on  the  similarity  variable  (t  u  cos  S/Lh). 

In  the  almost  sonic  case  (|cos  0|  «  5),  the  differential 
operator  becomes 

.3 


5  h      ,-  c_^  _  2L  ^ 
Sh^      ^t- 


(6.20)    D^f  1^  (u  ^  -  2L^) 


(at  this  point  y,  -x  could  replace  h,  t).   The  operator  in  the 
parenthesis  of  equation  (6.20)  does  not  seem  to  have  been 

studied.   We  can  expect  the  solution  to  depend  on  the  similarity 

3   /  2 
variable  t-^u  /h  L).  The  more  general  case  given  by  the  operator 

in  (6.18)  will  therefore  give  a  transition  from  one  similarity 
variable  (t  u^cos  o/lh) ,   to  the  other  (t-^u  /Lh  ) ,  as  we  move 
from  supersonic  flow  to  high  subsonic  flow. 

There  is  still  another  boundary  layer  to  be  considered. 
Looking  at  Figure  k.l,   we  see  a  singularity  H(y)6(x)  which 
corresponds  to  the  wake  behind  the  object.   The  three- 
dimensional  wake  has  already  been  considered  by  H.  Grad  [13]. 
To  consider  the  wake  we  have  to  "spread"  the  region  along 
the  positive  x-axis.   Carrying  on  the  same  analysis  as  used 
in  the  above  cases  vie   see  that  we  can  estimate  the  determinant 
D,  (6.2)  by 

(6.21)    D  ^  i  I  u^k^k^^  +  I  kg^L. 

1/2 
Here,  5  =  0(L   )  and  as  in  the  supersonic  Mach  layer,  the 

1/? 
thickness  of  the  wake  grows  as  (RL)    .   Corresponding  to 
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the  estimate  of  D,  (6.21)  we  have  the  differential  operator 

wl.dch  again  Is  essentially  the  diffusion  operator.   And  again 

we  can  expect  to  find  the  similarity  variable  (y  /hx)  to  occur 

In  the  solution  for  the  flow  field. 

In  our  discussion,  we  have  Ignored  the  case  when  the 

Mach  angle  is  small.   In  this  case  the  free  stream  velocity, 

u  ,  becomes  very  large.   This  is  the  hypersonic  limit,  and 
o 

we  do  not  consider  it  in  this  report. 


7.   Wake  Layer. 

In  this  section  we  apply  the  approximations  developed 
in  Section  6,  to  our  matrix  system  (5.4).   Recalling  that 
for  the  wake  we  had 

(7.1)  D  =  1  I  "oVa^  +  f  ^2^^ 
with 

(7.2)  5  =  0(L^/2) 

we  can  easily  compute  to  lowest  order  terms,  the  inverse 
matrix  of  system  (5"^)'   Performing  the  matrix  inversion  we 
get, 
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(7.3) 


l"^!^ 


_  ll    2    2 
"  D   "  3  ^2 


k 


fk^lL  u^k^k| 


u^k^-IlLk^ 


2     I^^o^i+IL^') 


liu^k^k^) 


2   -2 
-  3  '^'^o 


■^"o^l 


We  note  Immediately  from  the  matrix  in  (7.5)  that  the  last 
column  is  of  a  higher  order  than  the  second  column.   Therefore, 
the  solution  for  p*,  T,  and  V'^  is  to  lowest  order  terms,  only 
dependent  on  the  entropy  source.   To  this  order  then. 


(7.M 


P  =  -  ^ 


iu  T 
o 


5  k,-  ik^L/u^ 


and  on  Inverting  the  Fourier  transform,  the  density  is 
/„^\     ^        o      e     dk 


(27r)^j  k^^-  ik^L/u^ 


By  inspection  of  the  system  we  have  immediately 


(7.6) 


/J 


T  =  -p. 


so  that  the  perturbed  pressure  vanishes  to  this  order. 
Lastly  from  (7-3) 


(7.7) 


k-u  =  -^  iu  T 


5   °  k,-I.^„ 
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In  order  to  determine  the  velocity  field,  we  must  return  to 
the  original  system  of  equations  (5.4),  from  which  we  get 
to  lowest  order  terms 

(7.8)    (ivi-^^)(5)  =' (o°i -'£)''•"  ^ 


and  Inverting 
(7.9) 


'«!,   "^0°) 


r 


e^^'^  3k 


L4u  t 
o 


%\       (27r)2j  (k^-  Ik^L/u^)  '  5(27r)2 


(k^-  Ik^L/u)^ 


The  flow  field  Is  therefore  determined  by  the  Integrals 

1 


(7.10) 


and 


I  = 


r     e^^-^dk 


(27r)"  J    k^-   Ik^L/O; 


(7.11)  II  = 


1  r 


e^^-^       dk 


(27r)2  J    (k^-   lk^L/uj2 


"both  of  which  are  easily  evaluated. 

.    ,.        yx     \/2     -y^u  /4Lx 
(7.12)  i=M20(      O)  e 


2         '■ttLx' 


-2':' 


(7.13)        n  =  -2HM  (_^)^/^  e        ° 


Substituting  the  forms  I  and  II  Into  (7-5),  (7-6),  and  (7.9) 
we  determine  the  flow  field 
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o  ^  '^n  1/2  -y^u  /4Lx 

P  ^   1/2   -y^G  /4Lx 

(7.14)  ^         p. 

u   1/2   -y  u  /4Lx 

u„  =       7^  Tur  yH(x)^   '  v-inrfTrJ  e 


We  notice  from  the  solution  that  (y  /Lx)  Is  Indeed  the 
similarity  variable,  and  that  the  wake  thickness  is  0(£Lxy '^  ). 

At  this  point,  it  is  important  to  make  the  distinction 
between  solving  for  a  fundamental  solution,  and  solving 
for  flow  past  an  object.   In  the  former  case  it  is  not 
necessary  to  demand  that  products  of  the  solutions  are 
small  compared  with  respect  to  the  solutions.   On  the  other 
hand,  the  latter  case  absolutely  demands  that  this  be  true. 
This  is  in  order  to  agree  with  the  linearization  which  we 
performed  at  the  outset.   For  the  case  at  hand,  the  lineariza- 
tion  is  valid  when  the  quantity  tu^/^ttLx  is  small  compared 
with  unity. 

Regarding  (7.l4)  as  a  fundamental  solution,  and  taking 

the  limit  as  L  — >  0,  we  recover  the  fundamental  solutions  of 

the  Euler  approximation  (4.19-4.19).   This  follows  immediately 

by  noting  that 

u   1/2   -y2i^  /4LX 
(7.15)    lim   (th^)    e    °     =  5(y). 
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See  for  Instance  reference  [I?]- 

It  Is  Important  to  notice  that  the  convergence  to  the 
Dlrac  5 -function  is  not  uniform  with  respect  to  x.   For  any 
small  L  >  0,  we  can  find  a  large  enough  x  such  that  the 
wake  is  wide  enough  to  invalidate  an  inviscid  description. 
A  similar  effect  occurs  in  the  next  section  in  connection 
with  the  Mach  layer  spread.   We  will  further  elaborate  on 
this  point  there. 


8.   Mach  Layer. 

We  now  make  use  of  the  approximations  developed  in 
Section  6  for  the  layer  in  the  neighborhood  of  the  Mach  lines. 
Recalling  our  calculations  in  Section  6,  we  found  that  in 
order  to  obtain  a  description,  uniformly  valid  from  super- 
sonic to  high  subsonic  speeds,  we  should  take 

(8.1)     D  =  -  -I  iu  mn^sin  9  -  ^  ±u   m^n  cos   6   -  ^  Lm  . 
^  3   o  5   o  3 

1/2 
In  addition  it  was  found  that  6  varies  from  0(L  '  )  to 

0(L  '     )   as  the  flow  goes  from  supersonic  to  subsonic  speeds. 

Applying  this  to  the  Navier-Stokes  approximation  and  carrying 

only  lowest  orders,  we  find,  on  inverting  equation  (6,2), 

that 
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-m^sln^e        -m^ 


-urn  sin  6\    /    0 


(S.2) 


T 


1 

D 


Im2 


2  ^2 
-m 


-  -^u  m  sin  9 

;>  o 


1 ' ST  n    I  - 


u  m'sln 
o 


-u_m-^sin  9 


-|.»^ 


2  --2 

3  o 


u   rm   sin  0 
o 


Solving  for  the  transformed  density  to  the  lowest  order, 

,2.  . 


(S.JO 


P  = 


i(2u^+  5)t 
5u^L 


m 


2  "^ 

/n  sin  9        2mn  cos  0   21m-_N 

^   L        L      "  u.  '' 


and  then  inverting  to  get  the  density,  vje  have 


(8.4) 


P  = 


i(2u^4-  5)t 


^irnt  +  inh  ^^^  ^^ 


5  sin  0  L   (27r)  u 
On  inspection  we  see  that 


2                    5 
I  fU      .  2mn  cot  0   2im-'^csc  0\ 
o    (_  +  ^ ) 


u. 


o 


(8.5) 


so  that  to  this  order  the  flow  is  adiabatic.   Also  by 
inspection 


(8.6)     k-u  =  urn  Sim  9   p. 

As  before,  in  order  to  get  the  expression  for  the  velocity 
we  have  to  return  to  the  unseparated  system  (5-'^)>  ai^d  from 
there  we  get 


o 


/k. 


t(^q"J-  [(i+2Lu^.k)p  +  2Lk.u],j^;;-^ 


V 
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Again  only  carrying  lowest  order  terms  we  have 


(8.8)     (u,  X)  =  (-  i-,  m  cos  6  +  n   sin  0)  5  ^  , 
*"    '      ^1*2'^    ^   u^   u  m  sin  6  '    5  ^ 


On  taking  a  Fourier  Inversion  of  equation  (8.8)  we  get 

r 


(8.9) 


"l  = 


-1(2u^+5)t  ^ 


3u^sln  e  L  u^(2Tr)' 


y 


^Imt-Hlnh  ^  dm  dn 

— 2 5 

,r\      ,    2mn  cot   6     21iit'^csc   ds 

^L-  +  — i^ s: — ^ 


and 


r 


lu    (2u  +5)tcos   0  T 

(8.10)  {l2  =  -°       °  ^ 


5L 


^Imt   +  Inh  ^    ,^    ,„ 
e  m  dm  dn 


"o(2-)'J(£-% 


2mncot  0    21m''^csc  Q 


U- 


+ 


1(2G^+5)t  ^ 


3u^slneL  UQ(2Tr)2 


^ 


^imt   +  Inh  ^  ^  d^ 

2 ^ 

/n      ,    2mn  cot   0        21m  esc   0\ 
Vt~  +  f: 1? ) 


u 


o 


The  Integrals 
(8.11)    I  = 


u^(27r)2 


J 


Imt   +  Inh  ^  ,^  ,^ 
e         m  dm  dn 

— 2 3 

/n    2mn  cot  0   21m  esc  0  \ 


and 


/^ 


(8.12)     II  =  - 


u^(27r)' 


^Imt  +  Inh  „  ,^  ,„ 
e         n  dm  dn 
__ _ 

,n  2mn  cot  0  _  21m'^csc  0  s 

^L  "^    L  u^      '' 
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are  evaluated  in  Appendix  A,  in  closed  form  in  terms  of  an 
Airy  integral.   Therefore  with  the  integrals  I  and  II  evaluated, 
we  have  a  complete  description  of  the  structure  of  the  Mach 
layer  from  high  subsonic  through  supersonic  velocities. 
Inserting  these  forms,  the  uniform  solution  for  density  is, 

1  t\sin0  ^  4  t  ugcose   8  h  cos^g 

,^-2   ^N  ^'^        h^L  ^  L^sin   6  ^^  L  sin^0 

(2u  +5)t  e 

(8.13)      p'=-^ X 

hl/5  5  3inl/\2/3  3I/5  4I/5 


^  t2u^/53i^2g        ^1/6^^2/3,^3   e  h2/\2/342^^g 

^1(7^73  ^^273^73^273  -  L^/^V^sin^^^e  ^  ~P7§-T73J         ^^ 


where  A.  is  the  Airy  function  described  in  Appendix  A. 
Corresponding  forms  can  easily  be  written  for  the  remainder 
of  the  flow  field.  Because  of  the  complicated  forms  we  only 
exhibit  two  flow  regimes  which  are  expeclally  simple  to 
represent.   These  are  almost  sonic  flow,  when  |cos  9\    <   0(L), 
and  supersonic  flow,  when  both  cos  6   and  sin  0  are  bounded 
away  from  zero.   Note  though,  that  any  case  can  be  computed 
since  the  Airy  integrals  are  tabulated,  [11]. 

Supersonic  Flow 

In  Appendix  A  the  asymptotic  forms  of  the  integrals  I 
and  II  are  given  in  the  supersonic  limit,  (A. 25)  and  (A. 26) 
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respectively.  These  correspond  to  the  asymptotic  development 
of  our  Airy  integral  for  large  values  of  the  argument,  and  in 
our  present  notation  these  are  given  by 


(8.14) 
and 


1^ 


1   sin  e  l}^^   e 


-t  u  cos  9 
o 

hL 


4h 


1/2  „l/2  _.l/2«   1/^ 


u 


o 


cos  '      9    IT 


(8.15)    II  ^  0. 


The  latter  is  not  given  to  lowest  order  since  it  occurs  in 
combination  with  the  integral  I,  in  the  flow  field  representa- 
tion.  Inserting  these  calculations  into  equations  (8.4), 
(8.5),  (8.9),  (8.10),  we  have  to  lowest  order 

r 

"4lF 


t(2u^+  5)   u^cos  9     1/2  ■ 
P  = :: (— )    e 


10  u^cos  9 
o 


4Lh 


TT 


T  =  I  P 


(8.16)  <^ 


"i  = 


(2uJ+  5)t   u^cos  e   1/2 
'-^        /_0_ 


6u^cos  9 
o 


4TrhL 


) 


-t  u  cos  9 
o 

"Tlh 


-(2u^+  5)t   u  cos  9      1/2  ■ 

^   O        /  O        \ 

u^  =  — 7, ( )    e 


-t  u  cos  9 
o 

— 4lH 


v^ 


6ujsin  9 
o 


47rhL 


A  form  which  is  in  asymptotic  agreement  with  the  system  (8.16) 
already  exists  in  the  literatiire,  see  [l8],  [19]'   Examining 
the  flow  given  by  the  equations  given  in  (8.13)  we  see  that 
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the  remarks  made  in  Section  6  are  justified.   The  thickness 

1/2 
of  the  layer  varies  as  (hL)  '    ,   and  the  similarity  variable 

is  (t  u  cos  0/4Lh).   When  viewed  as  the  asymptotic  solution 

to  the  problem  of  gas  flcv;  past  an  object,  the  assiomed  linear- 

1/2 
ization  is  valid  when  T/(hL)  '   is  small  compared  to  unity. 

On  viewing  the  set  (8. 13)  as  a  fundamental  solution  and 

passing  to  the  limit  of  small  L  we  see  that  the  inviscid 

solutions  of  Section  4,  (4.17-4.19)  are  recovered.   It  is 

interesting  to  note  that  the  inviscid  solutions  are  also 

obtained  in  passing  to  the  limit  of  small  h.   This  limit  should 

be  thought  of  as  being  carried  out  for  a  fixed  but  very  small 

mean  free  path  L.   Physically,  this  process  tells  us  that  on 

approaching  the  boundary  layer  of  an  object,  we  approach  the 

conditions  of  inviscid  flow.   Again,  as  in  the  wake  of  Section 

6,  the  convergence  to  the  Dirac  6-function  is  not  uniform. 

For  any  small  L  >  0,  v;e  can  go  far  enough  away  from  the  origin 

such  that  the  flow  is  viscous.   As  a  criterion  for  viscous 

flow  we  can  choose  that  radial  distance  at  which  the  Mach 

layer  is  double  the  body  size.   A  simple  calculation  shows 

that  this  is  0(<3f  /l).   This  shows  that  even  though  the  flow  is 

Inviscid  outside  the  boundary  layer,  at  sufficiently  large 

distances  the  flow  field  again  must  be  described  by  a  viscous 

theory. 
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Sonic  Flow;        |cos0|  <  0(L) 

This  case,  as  mentioned,  also  allows  a  simple  represen- 
tation.  The  appendix  contains  the  representation  of  integrals 
I  and  II  in  the  sonic  limit.   In  this  limit  we  can  return  to 
the  x-y  coordinate  system.   From  (A. 20)  and  (A. 21)  we  get 


(8.17) 

and 
(8.18) 


3 

1  ^  ^c 


I  = 


IL^/^e  ''^'^ 


Al(- 


X  uzf  -^ 
o 


3I/3  4l/5|y|l/5  uj/5    3^^/^  2^/5  y^/5  l2/3 


II 


1  X 

3  y 


Equation  (8.18)  is  approximate  since  a  higher  order  term  is 
neglected.   The  latter  term  is  no  longer  negligible  for  large 
values  of  the  argument  of  the  Airy  integral  (outside  the  Mach 
layer).   If  a  study  of  the  full  flow  field  is  asked  for,  the 
remainder  term  must  be  included.   Inserting  the  above  forms 
into  (8.4,8.5,8.9),  we  get  the  sonic  flow  field. 


P 


{8.19)!  ^ 


/^ 


ri 


/ 


-1 


u. 


-u^x 

3y~y 


5t  e 


3 
1_^  ^ 

27  ?r 


x2u2/3 


V34I/5 1^,1/3^^/3^1/3  ^l(  4/322/34/3^2/3) 
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The  expansion  of  the  Airy  function  in  the  layer 
is  given  in  the  appendix;,  (A.  17).   In  the  above,  use  has 
been  made  of  the  relation,  u   =  5/3  which  holds  for  sonic 
flow.   Notice  that  u^  is  of  a  higher  order;  we  keep  it  in 
order  to  exhibit  it  in  its  lowest  order. 

Again  the  conclusions  drawn  in  Section  6  are  verified. 

2   l/^ 
The  thickness  of  the  layer  varies  as  (y  L)  ' -^   and  the  simi- 

larlty  variable  is  (x-^u  /y  L).   For  flow  past  an  object, 

T/(hL  )  '      must  be  small  compared  with  unity  in  order  to 

justify  the  assumed  linearization.   One  also  notices  that, 

irrespective  of  the  problem,  the  theory  Indicates  that  in 

order  for  the  linearization  to  be  correct,  one  must  be  further 

away  from  the  body  in  transonic  flow  than  in  supersonic  flow. 

In  order  to  examine  the  sonic  flow  field  in  the  limit  of 

vanishing  mean  free  path,  it  is  convenient  for  formal  reasons 

to  return  to  (8.4,  5,  9,  10).   From  these  relations  we  see 

that  the  sonic  flow  field  may  be  represented  by 


(8.20) 


5 


1 

^0 

Sx 

2 
5-0 

^x 

^x 
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v/here 

f      • 
^     ^-imx  +  my 

(8.21)    ^=-26 


(27r)^     („2_  ^2  im^j 


acts  as  the  velocity  potential.   Although  this  is  a  divergent 
integral  it  has  meaning  in  the  sense  of  distributions,  and 
the  evaluation  follows  that  of  Appendix  A.  This  yields 

(  .!!S 

?(  =  \y\{  —2 A,  ( 2 \   + 

1  5I/3  //3  ^2/5     '  3^/'  22/'  yV?  l2/5 
2^' 
3I/54I/3-I/3   -  ^        ^2^  2/3 

y2/3  l1/3  ^  5V3  22/3  y4/3  l2/5 

where  the  prime  denotes  differentiation  with  respect  to  the 
argument.   In  order  to  show  that  the  solution  (8.20)  becomes 
the  inviscid  sonic  solution  of  Section  4,  we  need  only  show 
that  the  bracket  of  equation  (8.22)  becomes  a  delta  function 
in  the  limit  of  zero  mean  free  path.  We  do  this  in  the  usual 
way,  by  letting  the  bracket  operatr  on  some  testing  fionction 
•i^(x).   Consider  the  Integral 

r 


(8.23)    I  ^^  i^(x)dx 

J    iyi 


taken  across  the  origin  between  some  definite  limits.  Under 
the  transformation 
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(8.24)    s  = 


^  1/5 
xu  '  ^ 
o 


^2/3  2I/3  L^/^  y2/3 


this  becomes 


r    2s- 


(8.25)    6 


se   -"  A^Cs'')^  ( 


1/3  ,2/3  pl/3  ^2/3 
^\.^:   (k ^—^ ^ s)ds  + 


+  6 


'"  -  I  s^  ,   2    .1/3  .2/3  pl/3  ^2/3 


J 


A^(s")t^  0 


u 


T/T 


Expanding  the  testing  function  in  a  Taylor  series  and  passing 
to  the  limit,  we  have  our  result.   Notice  that  again  the 
convergence  to  the  delta  function  is  not  uniform.   The  lowest 
order  term  after  the  delta  function  is  0(L  ^^   y  ^^).      As  in 
the  case  of  supersonic  flow  and  the  wake  before  that,  for 
any  small  but  fixed  L,   vje  can  go  far  enough  away  so  that  a 
viscous  theory  must  be  used.   Then  quadratic  terms  are  small 
compared  to  linear  terms  and  the  embarrassment  of  infinite 
drag  is  removed. 


9.   Contribution  of  Lift  and  Heat  Addition. 

In  Section  4  we  restricted  the  particle  source  so  that 
to  lowest  order  it  would  act  only  as  a  drag  singularity.   This 
excluded  the  possibility  of  lift  and  heat  addition.   In 
general  however,  these  two  have  to  be  included,  buth  for  the 
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asymptotic  flow  far  from  a  body  and  In  the  fundamental  solu- 
tion.  This  section  deals  with  the  effect  of  the  additional 
terms . 

Returning  to  the  source  integrals,  in  place  of  equation 
{hA)   we  define 


(9.1) 


and 


J 


e  d^  =  u 


Vl 


r  '^ 


(9.2)     J  ^^6  de  =  2qu/ 


Previously  the  heat  addition  and  lift  terms,  q  and  A,  v/here 
taken  as  zero.   Inserting  (9-1)  and  (9.2)  into  equation  (4.1) 


our  source  terms  become 
(9.3) 


2u' 
C^^*  — —  (q-T) 


2    3 

l,\\ 


(9.4) 


3  ■'. 


V 


'U. 


.2  /   " 


3  / 


A 


We  will  consider  the  effects  of  lift  and  heat  addition  in 
the  wake  layer  and  Mach  layer  separately. 

Wake ; 

The  source  term  in  the  contracted  system  (7.3)  now  has 
the  form 
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(9.5) 


2^ 


(q-x) 


k, TU  +  k^Au 
1   o    2  o 


From  inspection  of  equation  (7-3)  we  see  that  to  the  lowest 
order,  again  only  the  entropy  source  enters.   This  does  not 
mean  that  the  asymmetry  of  the  lift  is  Ineffectual  in  the 
wake.   As  we  will  see  shortly  the  lift  source  produces  an 
asymmetry  in  the  wake  velocity  field.   To  this  order  then 
the  density  field  is 


V       1/2        ..2 


(9.6) 


^  =  |^o('^-^)  (ir;;^) 


u 
■■JFttLx 


y^u^ALx 


and 
(9.7) 


T   =   -p. 


A  simple  calculation  shows  that  the  veolcity  field  is  given  by 
-1  =  -H(-)  (to)     ^ 


(9.8) 


V^ 


^   1/2   .y2-/4L^ 
-2  =  ^H(-)  (to)    ^    ° 

o   X   ^ttLx' 
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Mach  Layer; 

Returning  to  equation  (8.2)  we  have  to  replace  the 
source  term  by 


(9.9) 


V. 


-f  u^  (x-q) 

-u  m(T  sin  0  -  A  cos  6) 


A  simple  calculation  shows  that  with  (9'9)*  we  need  only 
replace  (2u  +  5)t  by 

(9.10)    |(2u^  +  5)t  -  2u^q  -  5>^   cot  ej 

in  all  the  formulae  to  get  the  additional  effect  of  lift  and 
heat  addition. 
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Appendix  A :  Evaluation  of  Integrals. 

In  Section  8,  there  occurred  the  two  related  Integrals, 


(A.l) 


I  = 


(27r)2u 


o  J 


^Imt  +  Inh  ^  ,^  ,^ 
e         m  dm  dn 

— 2 3 

/n  +  2Tnn  cot  9   _   21m  c_sc_9_\ 

^L       L  ul     ' 


(A. 2) 


II  = 


{2ir)\ 


^Imt  +  Inh  ^  ,^  ,„ 
e         n  dm  dn 
_ _ 

-n    2mn  cot  0   21m-^csc  9 


^L"  +  T 


u. 


Under  the  transformation. 


(A. 5)    I 


kuv 


m  = 


2u^H 


h  = 


t  = 


4u: 


2u, 


h' 


V 


(A.l)  and  (A. 2)  become,  respectively 


(A. 4) 


I  = 


(27r)' 


r    e^^^^'  ^^"^'  ^   dM.  dv 

~5 5 

V  +  M-v  cot  9   -   1^1,  CSC  9 


(A. 5) 


II  = 


27r' 


r  ^w  +  ivh'  ^  ^^  ^^ 

-^ 3 

V  +  iiv  cot  9   -   l|jL''^csc  9 


Integral  I. 

Since  h'  >  0,  we  can  complete  the  contour  In  the  upper 
half  complex  v-plane.   We  are  therefore  interested  In  those 
roots  of  the  denominator  which  have  a  positive  Imaginary 
part.   The  roots  of  the  denominator  are  given  by 
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(A. 6) 


V   +  = 


M.[cos   e  ±   iH[i   sin  9  +  cos^  ] 
2   sin  G 


Choosing  a  cut  just  below  the  negative  real  axis  of  the 
complex  |x-plane,  we  see  the  v  Is  always  In  the  upper  half 
plane  in  the  first  Rlemann  sheet.   The  Integration  in  the 
v-plane  gives 


(A. 7)     I  = 
Next  taking 


i  sin  e 
2-w 


(A. 8) 

we  get 
(A. 9) 


-co  Vcos^e  +  4i|x  sin  0 


2/5. 


\/cos^e  +  Vi  sin  0  =  ^  ^^^^y^  ^,  C 


I  = 


sin^^^e  6""^'^°^^^/^  ^^^  ^  +  h'cos^0/8  sin^e 


2Trh 


WT 


r 

I 

c 


-3_,;.2/t'sin  0   h'-'-'^^cos  0^    »  h'^'^^cos^0 
h''^'^^     2  sin^^^0       4  sin^'^^© 


X 


dC 


with  the  path  C  shown  in  Figure  A.l 


h*   ^cos  0 
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Finally  we  perform  the  translation, 

1/3 

,,  ,„s    .   s    1  /t'sin  0   h;_^^j^cos0N 

(A. 10)    C  =  -TT^  -  3  V 2/3 5/3-^ 

to  get 

2  t'^sine   ]^  t'^cos^e  2  h'cos^e 

.   .  2/3„   27   ui2   "9   sin  e   "^  27  „,„2o 
/.  ,,v   ,   i  sm  '-^0  ^  '    h'      ^  sin  0    .  /  \ 

A.(x)  is  an  Airy  integral,  given  by 
(A. 12)     A,(x)  =^    L^^/^--^dO, 

see  [11].   The  path  C  is  of  the  type  shown  in  Figure  A.l. 
For  the  case  at  hand 

/«  -,,\         1    /t'^sln  0     t'cos  0    _^  h!^_cos^x 

(A. 13)    X  =  -jr^^  ( jry^ r-A^ yTT'  ^   W/Tl — '' 

3V  ^   h'^/-'     h'-^^''sin-^/^0    sin^^^0 


Integral  II 

This  integration  follows  the  methods  used  in  the 
integration  of  I.   Carrying  out  the  integration  in  the 
complex  V -plane  we  get 

.       f^   itxt'    +  iv  h'  I P — 

(A. 14)  II  =  |-         e "^      (-COS0+YCOS   0   +  4i      sin  F )  d|a 

"°°  /cos^0  +  h\.\i  sin  0 


^      r^  int'    +  iv_j_h' 
2? 

'-  00 


=   -   cot  0  I  +  irr         e  "^     dM- 
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Considering  only  the  second  part, 

.   r°°  l|j.t'+iv  h' 
(A. 15)    II'  =  Ijj:  I   e       +  dn 

-00 

and  Introducing  the  changes  in  variables  given  by  (A. 8)  and 
(A. 10),  we  get 

1/3 
(^    i;c\    TT1    2       1       /t'slne   h'-''^  ^  cos  6      vt 
(A. 16)    II'  =-_ jw^ j-^  ( — ^^ ^^ — )I 

^  sln-^'^-'e  h'-^^^  h''^'^^      ^   sin'^'^^e 


2  t'^  ^.  fl  It'^^^^P  1  t^cos^e  ,  2  h'cos^e 

1/5   27  — 2  ^^"^"  Q  H"^  '^'^^  ^"  Q  sin  ^  "^^T  5 

21  sin-^^e^'^'  h'^       ^  "       y  sin  y   ^r  gin^e   „',  x 


Here,  x  is  again  given  by  (A. 15),  and  A.(x)  denotes  the 
derivative  of  the  Airy  integral,  (A. 12). 

Expansions  of  the  Airy  integral  may  be  found  in  [11] 
The  power  series  expansion  is  given  by 


3l/3('  2).  (X  +  |y  x^  . . . ) 


and  the  asymptotic  expansion  is  given  by 
(A.18)    A,{x)~^ -1^51^(1-1^1  +  ...) 

with 
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(A. 19)     ?  =|x^/2. 

In  connection  with  almost  sonic  flow,  the  integrals  may 

be  written  to  lowest  order  as 

2_  t'^ 

(^•20)   I  =  j-ri7T;T7*5  ^'^  ""'  ^(^) 

and 

(A. 21)    11^-  1^1, 

with 

(A. 22)    x=-l7^-^. 

Since  the  argument  of  A.  is  0(1)  the  power  series  expansion, 
(A. 17),  may  be  used  in  this  case.   In  connection  with  the 
supersonic  limit,  when  both  cos  9  and  sin  0  are  bounded  away 
from  zero,  we  consider  the  integrals  I  and  II  as  L  — >  0. 
Then  with 


(1^    o:^\          ^       h'^^^cos^e 
(A. 25)    X  -jn^ jr/T- 

3^^sin^/-^e 
and 

/„  oii^    £   2  h'cos-^e   1  t'cos^e  ,  5  t'^cose 
^'   sin^e    ^  sin  0    ^°    " 

1  t'^sine 

-TO  -;^ . 

we  get,  using  (A.I8),  that 
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(A. 25)  I^ 


717^ 


1  t'^cos  e 

1   sin  e  ^"  ^         FP 
1/2     1/2     ® 


2h'  cos    e  TT 


and 


(A. 26)  II  ^  0. 
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